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Abstract 

Modification of tlie particles in tlie course of the source evolution is consid- 
ered. Influence of this effect on multiplicities and correlations of the particles 
is displayed, including an enhancement of the production rates and identical 
particle correlations and also back-to-back particle-antiparticle correlations. 



1 Introduction 



It is well known that in quantum cliromodynamics the chiral invariance is 
spontaneously broken at low energies and the fermion masses (say nucleon or 
constituent quark mass) arise essentially due to chiral simmetry breaking. In 
the course of phase transition at high fermion density or at high temperature 
the approximate chiral invariance is restored and the light {u,d,s) quark 
masses become small (ensuring small explicit violation of the symmetry). 
In practice the change of the fermion masses is treated (semi) classically. 
However it appears that, if the change of the mass is fast enough, then the 
quantum fermion field is changing in a specific way producing correlated 
fermion-antifermion pairs. 

Similar effects were considered earlier for bosons (mesons [0, 0, ^ and 
photons p) with application to heavy- ion colhsions p|, |[. Below we consider 
the case of spatially homogeneous large volume. Applications of this effect 
to heavy ion collisions will be given elsewhere. 

Decomposition of the free fermion field is taken in the form: 

V^(^) = j (^e^''' EK(k)6.(k)e-^* + ^;.(-k)4(-k)e^^*] (1) 

where hu,h\, and d^,dl are annihilation and creation operators of particles 
and antiparticles, obeying standard anticommutation relations. 

Bispinors u^(k),v^(k) are orthonormal, and bispinors u^(k) and f^(— k) 
having opposite momenta are orthogonal. For these bispinors related to Dirac 
equation 

(^7°^-7"A;'^-mj^(k,t)=0, n = l,2,3 (2) 
we take the standard representation 

Mk) = ^(m + 7A:)(Q), ^^(k) = ^(^-7^)(°) (3) 

with two-component unit spinors {si)p = 6ip, {s2)p = and normalization 
factor 

N = N{E,m) = {2E{E + m))^^^ (4) 
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2 Step-like variation of the fermion mass 



Let the mass m in (5) depends on time, m = m{t). We consider in this 
section the step-hke variation of the mass. This simple case reveals the main 
features of the phenomenon. Let at time t — Q the mass changes instantly 
from rrii to m/. According to Dirac equation (2), the time derivative of 
the field ip has jump discontinuity at t = 0. Therefore the field '^(k, t) is 
continuous at the point t = , 

iji{\i,-5t)^ijf{\i,+5t), St^ (5) 

Turning to decomposition (1) we choose the third coordinate axis along the 
direction of momentum k: 

A;(3) = A; = |k| (6) 

and use for u^,v„ the representation (3). In this case the bispinors u„,Vi, 
correspond to definite helicities /i (/i = 1/2 for i/ = 1 and h = —1/2 for 
v — 2) and Eq.(l) gives decomposition over states with definite momentum 
and helicity. In view of continuity of the field ip we have 

Mi,^(k)6i,^(k) + Vi^^{-k)dl^{-k) = Uf,MbfA^) + VfA-^)d\A-^) (^) 

(sum over z/). Multiplying both sides of (7) by ■uj-^(k) and then by fJ-^(— k) 
and using orthogonality conditions we get the Bogoliubov transformation 
of the annihilation and creation operators expressing final-state operators 
through initial-state operators: 

6/,,(k) = a,(k)6,,,(k) +/3,(k)4,(-k), 
4,(-k) = -/?,(k)6,,,(k) + a,(k)<,(-k) (8) 

for every u, k, with coefficients 

a„(k) = 4,„(k).,„(k) = 4„(-k),v(-k) = (E, + m,m-,,n,) + e ^ 
Mk) = t.},„(k).,,„(-k) = -4„(-k)«,,„(k) = ^ m + m,-E,-m>) 

^ / 

for z/ = 1, 2 correspondingly where A^'^, Nj are normalization factors of the 
spinors u,y, in the initial and final states, given by (4). As can be seen from 
(9), the Bogoliubov coefficients a,y(k), /5j,(k) satisfy the condition: 

|a4k)|2 + |/3.(k)|2 = l (10) 
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so that the transformation (8) is unitary SU(2)-transformation, as it must 
be to preserve the anticommutation relations. In the case (6) the Bogohubov 
transformation connects creation and annihilation operators which have op- 
posite directions of momenta and opposite spin projections (equal helicities). 
It can be shown that the last statement remains valid for any direction of the 
momentum k. Let us note that for step-like transition the Bogoliubov coeffi- 
cients (9) are real- valued. The coefficient /5,y(k) changes its sign if we change 
the sign of the momentum (k — k), or change 2), or interchange 

the initial and the final states [f ^ i). 

In general the coefficients a^(k),P^(k) of 5'C/(2)-transformation (8) can 
be represented in the form 

Q;(k) = cosr(k)e^^", (3(k) = smr{k)e''^^ (11) 

where 

r{k) ^tan-^\p/a\ (12) 

is the main evolution parameter and the phases ipa, do not play important 
role and they will not be considered here. For step-like transition the phases 
are absent. 

Using (6), (9) we get 

As one can see from (13) the evolution parameter r(k) is equal zero at k = 0, 
it is maximal at 

2 2 

e ^kl^mimf+ . "^'""^ (14) 
{rrii-mfY 

and it falls down slowly at large A;, 

\mi-mf\ 

^(k) ~ (15) 

In the point of the maximum the parameter r depends only on the ratio 
mf/rrii and it is sizable if the mass ratio is small (or large), say ^ m/. 
In this case 



tanr = |/3/q;| 1 — 2w— for k'^^rriimf (16) 
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If one takes the fermion masses rrii and mj to be the constituent quark 
mass(~ 350 MeV) and the current quark mass (~ 5 MeV) then tanr^Qj. ^ 
0.76 at ^ 42 MeV. 

In reality the change of the fermion mass has finite time duration r and 
the step-like approximation is valid only if the momentum is much less than 
inverse time duration, A;r ^ 1. The effect of finite time duration will be 
discussed in the next section. 



3 Smooth transition 

Let us consider the smooth variation of the fermion mass m{t). In this case 
the coefficients a(k), /?(k) of the Bogoliubov transformation can be expressed 
through solutions of the Dirac equation (2) ||]. To solve the equation it is 
helpful to use its squared form, representing the Dirac ^/'-function in (5) in 
the form 

ij{k, t) = (^^7°| -lk + m(t) j x(k, t) (17) 
Then the equation (2) takes the form: 

(^^,-^l'^ + k^ + m\t)'jx{k,t) = (18) 

which splits into two complex-conjugated equations for two upper and two 
lower components of x- The bispinor x can be written in the form 

X.(k,t)= ('^^)^(k,t)+(°J^(k,t), z. = l,2 (19) 

where are unit two-component spinors, see (3). As a result we have the 
second-order equation for scalar function ip{]<.,t): 

which is nothing but oscillator equation with complex-valued variable fre- 
quency (energy). 

Considering the smooth variation of the mass m(t) we use the reference 
model with 

m{t) = ^^^l±^ + ("^^p-) tanhi2t/r) (21) 
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for which (20) has the exact solution and contains the important parameter 
T giving characteristic time of the mass variation. The resulting solution (17) 
of the classical Dirac equation (2) describes the evolution of the normalized 
fermion field 0. Choosing one of the upper spinors (say si) we get 

^{k, t) = Mi(k, m„ Ei)e-''^'\ t < 
V'(k, t) = ai(k)ui(k, mj, Ef)e-'^f^ + /3i(k)ui(k, mj, -Ef)e''^f\ t > (22) 

showing appearance of the negative-energy wave (creation of antiparticles) 
with correct Boboliubov coefficient The choice of another unit spinor 
in (22) gives similar result. The momentum k in (22) may have arbitrary 
direction. If we take the k-direction along the third (spin quantization) axis 
then bispinors Ui in (22) have definite helicity h = 1/2 which is conserved in 
the course of transition. For the model (21) we get: 



ll/3(k, 



sinh [^{±Ef + Ei-mf + rrii) 


) sinh (^^{Ef±Ei±mfT mi)) 


sinh {^Ej) 


sinh (f E,) 
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The Bogoliubov coefficients in (23) satisfy the crucial unitarity condition 
(10). In the limit of step-like transition they correspond to (9) as one can 
check by straightforward calculation. 

As can be seen from (23), the coefficient /?(k, r) — ^ at k — and it 
falls down for large k: 

1^/, M /, N sinh (Trrlmj — m/-|/4) , . ^. 

In the region rrif <^ k <^ rrii, where the coefficient /3(k) and the evolution 
parameter r(k) are maximal for step-like transition, we get: 

l^(k..)P«i-i(^ + ;|M™.r/2)) (25) 

with 

h{x) = -{1 + xcothx), h>l (26) 

where rrii ^ mf,kT ^ 1 was taken in accordance with typical parameters 
rrii ~ 350MeV, rrif 5MeV, k ~ AOMeV in Section 2 and with typical time 



(23) 
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duration r ~ 1/m. As a result the maximum of /3(k, r) and r(k, r) is reduced 
for finite r and shifted to smaller momentum, 



|/?max(k, r)|^ ^ 2 ~ \/^"^//"^«' ~ ^JmiUif/h (27) 

due to /i-factor. At large A; ^ r^^ the effect is exponentially small, \(3\ ~ 
exp(— 7rA;r/2) for smooth transition. 

4 Fermion creation and their correlations 

Using Bogoliubov transformation (17) one can find the final-state fermion 
momentum distributions and the final-state fermion correlations. We confine 
ourselves to symmetric case when the fermions with opposite momenta are 
produced in an equivalent way (say central collisions of identical nuclei). For 
simplicity the Bogoliubov coefficients will be taken to be real-valued and 
= |k| dependent. We also consider the simple model - fast simultaneous 
transition of large homogeneous system at rest. The movement of the system 
can be taken into account by shifting of each moving element to its rest frame 
and integrating over proper times and space-time raidities of the elements of 
the system as it was done for photon production in heavy ion collisions. 

The final-state momentum distribution of the fermions (single-particle in- 
clusive cross-sections) is given by (below we use notations hy, dy for operators 
having definite helicity = ±1/2) 

^^/,.(k) = Y^ = (4,.(k)&/,.(k)) (28) 
where brackets mean averaging over initial state. Using (8) we get: 

Nj^M = |a(k)p(6t,(k)6,,.(k)) - \(3{k)\'{dU-mA-^)) + l/^^l'^^ 

(29) 

and similar expression for antifermions (with interchange h ^ d). The last 
term in rhs of (29) reflects the result of the vacuum rearrangement (the 
initial vacuum of fermions having mass mj is not the ground state of bf,df- 
operators). The factor V/{27t)^ replaces 5^(0) in the case of large but finite 
volume V. 
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It is convenient to introduce the level population function n,^(k): 

^.(k) = J^sM^) (30) 

Then (29) and corresponding equation for antiparticles take the simple form 
(independently for every helicity) 

n/,.(k) = |a(k)|V.(k) + |/3(k)|2(l -ri,,.(-k)) , (31) 

%,.(k) = |a(k)pn,,.(k) + |/3(k)|2 (1 - n,,.(-k)) (32) 

where the notation n is used for antifermions. In general one has to consider 
simultaneously the fermions and antifermions with momenta ±k to get the 
full picture of the particle creation. For example, if there are no particles in 
the initial state then wc have |/?(k)p particles of each kind in the final state. 
If there is one particle having momentum k and helicity in the initial state 
then we have 

n/,,(k) = l, %,,(-k) = (33) 

and particles in the rest of the states. If there is a fermion-antifermion 
(singlet) pair having zero momentum in the initial state, that is nj^,^(k) = 
1) ^1,1/ (~k) = 1, then we get 

nf,^{k) = n^,.(-k) = |a(k)|^ n^,.(-k) = n^,.(k) = (34) 

for the final state. Therefore in zero-momentum states which are initially 
copletely occupied (both with fermions and antifermions) the occupation 
number decreases contrary to initially empty states where the occupation 
number increases. For termal equilibrium n(k) is the Fermi distribution, 
depending on temperature and chemical potential and the both effects are 
present. 

The transition effect is better seen in two-particle correlations. Two- 
particle inclusive cross-sections are given by 

{6}^„(k2)6,,„(k2))(6t,^(k06,,„(k0)-i„„{t},„(k2)();,„(k,))((.},^(k,)(,/,„(k2)>(35) 
for two fermions with helicitics /i, u (correlation of identical fermions) and 

(6t,(k2)6^,.(k2))(4^(k0c^y,^(k0) + V(4,.(k2)4,(ki))(ci^,^(k06^,.(k2))(36) 
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for fermion-antifermion pair. The last term in rhs of (36) is essential only in 
the presence of the time evolution of the fermion field. Here for simplicity we 
do not take into account the dynamical correlations arising due to Coulomb 
and strong interactions of the fermions. 

Considering correlations in finite volume V we can use the modified cre- 
ation and annihilation operators ^ ^ which are nonzero only in the volume 
V and satisfy modified anticommutation relations of the form 

6^(ki),6t(k2)]^ = V(^G'(ki-k2), G(0) = 1 (37) 

where G(k) is the normalized form-factor of the fermion source (normalized 
Fourier transform of the source density). The corresponding correlator of the 
modified creation and annihilation operators has the form 

(6t (ki), K{k2)) ^ 5,.^n((ki + k2)/2)G(ki - k^) (38) 

The right hand side of (37), (38) represents the smeared 5-function having 
the width of the order of the inverse size of the source which is much less 
than the width of the momentum distribution n(k). 

Using Bogoliubov transformation (8), substituting (38) for arising fermion 
correlators of initial-state operators and taking into account the small width 
of the form-factor G'(k) one can express the final-state correlators through 
Bogoliubov coefficients a,f3 and the form-factor G: 

= ^ (|«(k)pn,.(k) + |/3(k)p(l - n,,.(k)))' \G{k, - k,)p, (39) 

{blMdUki)){dfAki)bfAk2)) 
^ 'a(k)/5(k)(l-ni,,(k)-n,,,(k))G'(ki + k2)|2 (40) 



{2tt 

where smooth functions a(k),/?(k) can be evaluated at any of momenta 
ki, k2 ~ ±k and we suggest that the process is k ^ — k symmetric. Taking 
sum over helicities u we finally obtain the relative (divided by the product 
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of single-particle distributions) correlation functions which are measured in 
experiment: 

c++(kl,k2) = l-^lG(kl-k2)|^ (41) 

C+-(ki, k2) = 1 + ii?2(k)|G(ki + k2)|2 (42) 

with 





|a(k)/?(k)(l-n,(k)- 




2 


[|a(k)Mk) + |/3(k)|2(l-n,(k))] [|a(k)| 


|%(k) + |/3(k)|2(l-n,(k))] 



(43) 

where we suggest that the helicities ±1/2 are equally represented in the initial 
state. 

The equation (41) represents the identical particle correlation in its sim- 
plest (no interaction) form. The minus sign is a distinctive feature of the 
fermions and the factor 1/2 arises due to sum over polarizations (only two 
of four final polarization states contain identical particles). 

The equation (42) describes the time evolution effect which depends on 
evolution parameter r (or (3, a). For the vacuum initial state the relative 
correlation function is big being equal to c? 120^ > 1/2. It is especially big 
for small f3 but in this case the pair production itself is weak, see (31), (32). 
In the case of cold dense matter the correlation function is big for empty 
(high momentum) initial states. In the case of hot quark-gluon plasma in 
the initial state one has to take into account the presence of effective termal 
quark mass which is much bigger than original quark mass (mj ~ gT for 
light quarks where g is QCD coupling constant and T is the temperature). 

5 Conclusion 

Calculation of the evolution effect for fermions shows that opposite-side 
fermion-antifermion correlations can be large. They can serve as a sign of 
chiral phase transition in quantum chromodynamics. 

The work was supported by the Russian Fund for Basic Research (grant 
02-02-16779). 
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